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ABSTRACT 

The Advanced Amphibious Assault Vehicle (AAAV) is a high water speed 
amphibious armored personnel carrier that will replace the current family of Marine 
Corps amphibious assault vehicles. The AAAV is currently in Phase I of the DOD 
Acquisition Process. During this phase-extensive development, testing 1s being 
conducted and prototypes are being constructed. Ongoing tests of the current electrical 
system architecture are revealing problems and issues that need to be addressed. Present 
testing is also revealing the need for a troubleshooting tool that can be used to simulate 
and test proposed solutions. 

A distributive computing effort is being conducted with Purdue University in 
order to provide the Marine Corps a tool where by it can test and evaluate the entire 
AAAV electrical system architecture. This document provides a general system 
description of the AAAV, an electrical system architecture overview, and a candidate 
electrical system description. Testing, modeling, and computer simulation efforts applied 
to the NBC (Nuclear, Biological, and Chemical) fan/filter motor drive and battery banks 
are described and the resulting data analyzed. Current research is investigating an issue 
regarding energy being forced back into the high power distribution bus due to 
regeneration from sudden impacts on the vehicle turret. An ultra capacitor bank may be 
incorporated in the AAAV electrical system in an effort to mitigate the effect of this 
regenerative energy. Therefore, testing, modeling, and computer simulation of an ultra 


capacitor is also reported. 
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I. INTRODUCTION 


A. AAAV Program 


The AAAV (Advanced Amphibious Assault Vehicle) will provide the Marine 
Corps a weapons system fully capable of implementing ship-to-objective maneuver as an 
integral part of the amphibious triad (AAAV, V-22, LCAC) to execute the concepts of 
Operational Maneuver From The Sea (OMFTS) and Ship to Objective Maneuver 
(STOM). Battlespace dominance by Marine forces will be significantly enhanced asa 
result of the AAAV’s high water speed and superior land mobility which have 
historically limited the rapid maneuver of armored combat vehicles. The AAAV will 
travel as fast as 29 miles-per-hour (mph) on water as opposed to the 10 mph the current 
AAV can travel. The AAAV is designed to allow immediate, high-speed maneuver of 
Marine infantry units as they emerge from attack positions aboard ships located beyond 
the visual horizon -- 25 miles and beyond. Projection of these forces will be conducted 
as a Single, seamless stroke that capitalizes on the intervening sea and land terrain to 
achieve surprise and rapidly exploit weak points in enemy littoral defenses. In Т" Ше 
AAAV is faster, can travel greater distances before refueling, 15 navigationally smarter 
and more maneuverable than its predecessor. 

The AAAV is currently in Phase I (Program Definition and Risk Reduction) of 
the Department of Defense (DoD) Acquisition Process. Phase II (Engineering and 
Manufacturing Development) is scheduled to begin Fiscal Year 2001 (FY01). The final 


phase, Phase III (Production, Fielding/Deployment and Organizational Support), is 


scheduled to begin FY06. The total cost of the program is expected to be $7.5B with 


each unit costing $6.6M. The Marine Corps is scheduled to receive 1013 units. [Ref.1] 


B. Research Focus 


A distributive computing effort is being conducted with Purdue University in 
order to provide the Marine Corps a tool where by it can test and evaluate the entire 
AAAV electrical system architecture. Reporting in this thesis includes a general system 
description of the AAAV, an electrical system architecture overview, and a candidate 
electrical system description. Efforts are then applied to documenting the testing, 
modeling, and computer simulation of the NBC (Nuclear, Biological, and Chemical) 
fan/filter motor drive and battery banks. Furthermore, research will also address the 
requirement to capture energy being forced back into the high power distribution unit due 
to regeneration from sudden impacts on the vehicle turret. An ultra capacitor bank may 
be incorporated in the AAAV electrical system in an effort to capture regenerative 
energy. Therefore, testing, modeling, and computer simulation of an ultra capacitor is 


also reported. 


С: Chapter Overview 


Chapter I of this thesis provides an introduction to the AAAV program. It also 
identifies the research focus and provides a chapter overview. Chapter II contains a 
comprehensive system description of the AAAV. It details a general system description 
and introduces the four modes of operation. In addition the discussion within Chapter II 


addresses the electrical system architecture, beginning with a topology overview and 


ending with a primary components description. Chapter III contains an introduction to 
distributed computing and background information on the programming language ACSL 
(Advanced Continuous Simulation Language), which will be used in this thesis to 
develop an NBC fan drive model. Further, the Automated State Model Generation 
Algorithm used by P. C. Krause & Associates (PCK&A), and Purdue University, to 
develop their NBC fan drive model is described and referenced in Chapter III. Chapter 
IV includes a description of the selected candidate electrical system architecture. In 
addition it provides background on the major components of the selected architecture, 
specifically the NBC fan, ultra capacitor, and batteries. Chapter V provides discussion 
regarding the modeling efforts applied to the NBC fan drive, ultra capacitor, and battery 
unit. It also provides results and conclusions for each model. Chapter VI is the 
conclusion. It summarizes the research while documenting results of the work. It also 


identifies Department of Defense relevance of the work and ideas for future work. 


THIS PAGE INTENTIONALLY LEFT BLANK 


П. AAAV SYSTEM DESCRIPTION 


A. Introduction 


This chapter provides a comprehensive system description of the AAAV. It 
details a general system description and introduces the four modes of operation. In 
addition the chapter addresses the electrical system architecture, beginning with a 


topology overview and ending with a brief description of the primary components. 


B. System Description 


The Advanced Amphibious Assault Vehicle (AAAV) is a high water speed 
amphibious armored personnel carrier that will replace the current family of Marine 
Corps amphibious assault vehicles, the AAV7A1 series. Operationally configured, the 
AAAV will weigh about 37 tons, be able to carry 17 combat-equipped Marines and a 
crew of three over 3 foot high waves in excess of 20 knots, and travel over land at 45 
mph. Armed with a 7.62 mm machine gun and a 30 mm cannon, the AAAV will use 
Global Positioning System, forward-looking infrared, and a night vision system for 
navigation, targeting, and intelligence gathering. [Ref. 1] 

The AAAV force is designed to provide a tactical assault and sustainment 
capability for Marines of the Marine Air-Ground Task Force. From amphibious ships 
standing well offshore, even over the horizon from the objective, the AAAV rapidly 
transports the landing force over the beachhead to an objective ashore, using maneuver 
and speed, plus on-board firepower to achieve superiority over enemy forces. Once 


ashore, the AAAV will serve as an armored personnel carrier, providing transportation, 


protection, direct fire support and command, control and communications (C3I) for the 
Marines. 

High water speed requires minimal hydrodynamic drag; to accomplish this, the 
AAAV will retract its tracks and suspension system and deploy cover plates on its 
underside to present a smooth surface. Thus, while in the water and traveling at high 
speed, the AAAV will be a planing hull craft. When the AAAV approaches 15 to 20 feet 
of water depth at the shore, it will reduce speed to come off plane, retract the cover 
plates, lower the suspension and tracks, and travel the remaining distance to the shore at a 
speed of approximately 8 to 10 knots. 

Once ashore, the Marine Corps will use the AAAV (P) as a tracked armored 
personnel carrier to accomplish dominant maneuver and precision engagement; in that 
mode, the 17 infantrymen carried inside will typically dismount to fight. The crew will 
use the AAAV (P)'s primary and secondary weapon (predicted to be a 25-35 mm cannon 
and 7.62 mm coaxial mounted machine-gun, respectively, with fully stabilized turret and 
forward looking infrared optics) to support the infantry and armor combat elements. 

The AAAV's land mobility characteristics must be comparable to the Marine 
Corps’s M1A1 Abrams main battle tank. This requires a top speed of approximately 45 
mph, the capability to traverse the same terrain as the tank during cross-country 
operations, and the capability to cross the same obstacles and terrain features (for 


example, trenches, hills, walls, and soft soils) as the tank. 
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Figure 2-1, Advanced Amphibious Assault Vehicle. [Ref. 1] 


System Operating Modes 


a. High Speed Water Mode 


An exploded view ofthe AAAV in high speed water mode is 
provided in Figure 2-2. The AAAV will operate in this mode 
during waterborne operations and will possess the following 


attributes. 


e Top Speed: 20-25 knots in 3 feet significant wave height 
e Maintains a heading in 10 feet significant wave height 

e Self-righting (100 degrees) 

e High water speed turning diameter <170 meters 

e Stops forward motion within 75 meters 


e Water range: 65 nautical miles 


Extended Bow Flap 


Retracted Suspension & Track 







Power Transfer Module 


Figure 2-2, AAAV High Speed Water Mode. [Ref. 1] 


Transition Mode 


A view of the vehicle in transition mode is documented in Figure 
2-3. The AAAV will operate in this mode while entering and 


exiting the water (from ship to land). Additional details include 


e Used to enter and exit the water (from ship to land) 
e Slow waterspeed operations 

- Top speed: 9 knots (4.5 knots reverse) 

- Surf transit: 8-foot plunging surf 

- Survivability: 360° rollover, 180° pitch-poling 


- Fully submersible for 5 seconds 


e Also used for riverine and coastal operations 








Power Transfer Module 


. 23" Waterjets 


Transition Flap 
Extended 


Suspension & Track HU 
Deployed and Powered 


Forward & Reverse 
Steering 


Figure 2-3, AAAV Transition Mode. [Ref. 1] 


с. Land Mode 


The vehicle in land mode is pictured in Figure 2-4. The AAAV 
will operate in this mode on land and possess the following 


attributes. 


e Cross country mobility: > = M1A1 Abrams tank 

e Top speed: 45 mph (0-20 mph in 7 seconds) 

« Range: 400 miles (land only) or 250 miles after 1-hour (25 
nautical miles) high water speed transit 

e Obstacle crossing 


-  8-foot trench span, 3-foot vertical wall 


Transmission 
:6 Speed 







Support Rollers 


Automatic 
Track 
Tensioners 


Suspension Units 


Double-Pin Aluminum Tra 


Figure 2-4, AAAV Land Mode. [Ref. 1] 


d. Silent Watch Mode 
A number of the key system components associated with the silent 
watch mode are illustrated in Figure 2-5. The AAAV will operate 


in this mode while stationary on land. In this mode, the following 


attributes hold. 


e Used for stationary vehicle to reduce acoustic and thermal 


Signature 
e Fully combat ready 


е One or more crewman on watch 
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- All com/nav/fire control/vision subsystems fully 


operational 


e Achieve full land mode capability within 60 seconds 


| m Whip Antennas 
„ ll [ 


i 





Conformal Antennas mw 


Weapons Station & Sight Systems PLE APU 
= Night Vision 2 ET. Diesel 
i « 330 Amps 


x Communications 
* SIRGARS SIP 21 (VHF) 
* SINGARS SIP 22 (VHF) 
= SINGARS ICOM (VHF) 
* AN/PSC-5 (VHF, UHF LOS, UHF SATCOM) 
* EPLRS 
« Wireless HeadsetAntercom 


- Batteries 
* 4 hour of operation 


Figure 2-5, AAAV Silent Watch Mode. [Ref. 1] 


C. Electrical System Architecture 


1: Topology Overview 
A block diagram of the electrical system 1s shown in Figure 2-6. This architecture 
schematic 1s of a level of detail greater than will be initially modeled by P.C. Krause & 


Associates or addressed in this thesis. The candidate electrical system for modeling and 


simulation can be found in Chapter IV, Figure 4-1. While a detailed description of the 
entire electrical architecture 1s outside the scope of this thesis, the primary components 


are identified below. 


AAAV Power Bus Architecture 


Vetranics 
Batteries 


Power 
Conholler 





=m + 


Figure 2-6, Electrical System Architecture. [Ref. 2] 


Primary Components 


a. Batteries 


- Eight (8), 12 volt, manufactured by Optima 
e Three series pairs dedicated for starting. 


e One series pair dedicated to the Vetronics Bus during starting 


b. Generators 


e Main generator — 800 Amp, six-phase wound-rotor 
synchronous generator connected to two six-pulse diode bridge 
rectifiers whose dc outputs are connected in parallel. The field 
windings are connected to an exciter that is a three-phase 
machine whose armature voltages are rectified using a rotating 
six-pulse diode bridge rectifier. 

e Auxiliary Generator — 330 Amp, connected to the 


Vetronics Bus 


C Ultra Capacitor 
e 100 Volt, 1F, 5000J capacitor manufactured by Pinnacle 
Research Institute. Not shown in Figure 2-6. See Figure 4-1. 
d. Hull Power Distribution Unit (HPDU) 


e Provides main power distribution to the Vetronics bus and high 


power bus. 
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Remote Acquisition and Control Modules (RACMs) 


• Five (5-HRACM) hull and one (1- TRACM) turret 
e Provide Secondary power distribution, Sensor Data Collection, 


Sensor Power, and Communications through Utility Data Bus 


Sealed Slip Rings 


• Provides power and data signal path to electrical systems 


located in the turret. 


NBC Fan Drive 


e 3.5 Hp, Wye wound, permanent magnet, brushless DC motor 
manufactured by Invincible Airflow Systems. Powered using a 


120-degree voltage source inverter with sensorless control. 


Ш. COMPUTER SIMULATION 


A. Introduction 


Chapter III contains an introduction to the theory behind the application of 
distributed computing to power system simulations. It also provides an overview of the 
programming language ACSL (Advanced Continuous Simulation Language), which will 
be used in this thesis to develop an NBC fan drive model. Further, it will introduce the 
Automated State Model Generation Algorithm used by PCK&A, and Purdue University, 


to develop their NBC fan drive model. 


B. Distributed Computing 


The consulting firm, P.C. Krause and Associates (PCK&A), has developed a 
technique that will aid in the simulation of complex systems in general and power 
electronic based systems in particular. This technique has been termed distributed 
computing. Distributed computing is a heterogeneous simulation environment that 
makes it possible to divide the simulation of power electronic based systems into any 
number of synchronized personal computers (PC). Moreover, subsystems can be 
implemented using different languages (1.e., Simulink, Saber, and ACSL) making it 
possible to have the advantage of a very large computer by paralleling inexpensive PC's. 

In addition to the gain in computational speed, another advantage of this 
simulation environment is the ability to decouple the user interface and the simulation 
that 1s controlled by the interface. In particular, the user interface can be executed on any 


computer that is connected to the Internet allowing users at remote locations to “exercise” 


the simulation and retrieve data. Also vendors who develop proprietary simulations can 


link the model as executable and not divulge the model inner workings. 


C. ACSL Background 


The acronym ACSL (pronounced “axel”) stands for the Advanced Continuous 
Simulation Language. It has been developed expressly for the purpose of modeling 
dynamic systems described by time dependent, nonlinear differential equations, and/or 
transfer functions [Ref. 3]. Simulation of physical systems is a standard analysis tool 
used in the evaluation of hardware design prior to and in concert with actual construction. 
An accurate simulation model may provide invaluable insight into design flaws or system 
limitations before extensive capital investment is made on hardware. For instance in 
analyzing the AAAV electrical power system, the engineers of General Dynamics 
uncovered an instability in the electrical system during low power operation when the 
batteries are disconnected. Discovering the instability early will allow the General 
Dynamics team to troubleshoot the instability. Knowing this up front, the proposed 
design could then be modified without costly retrofits or adhoc solutions. 

For the power engineer, ACSL facilitates the analysis and control of complex 
electromechanical networks. Variable speed electric drives are popularly used in 
pumping applications (fans, blowers, and compressors), large machinery actuation 
(elevators, hoists, cranes, conveyors, and crushers), remote servoing (robotics and 
precision machine tools), electric vehicle propulsion (cars, busses, trains, and ships), and 
as an alternative to hydraulic and pneumatic actuation systems. In addition to the 


electromechanical device, all of these systems are characterized by nonlinearities which 
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make application of linear system theory difficult (if not impossible). After designing 
comprehensive electric drives containing power electronic converters, electric machines, 
multi-loop analog and digital controls, and various protection and sensing devices, it is 
advantageous for the engineer to develop computer-based models which predict as 
accurately as possible the operational modes of the interconnected system. As the 
hardware becomes available, the models are validated and refined as necessary. Thus, 
modeling and simulation are integral parts of the design and analysis process. 

The ACSL program is intended to provide a simple method of representing these 
mathematical models on a digital computer. Working from an equation description of the 
problem (with logical conditions conveniently handled), the user writes ACSL statements 
to describe the system under investigation. Run-time commands exercise the model; 
these statements can be submitted for solution in command files or entered interactively. 
Simulation studies may be conducted, variables plotted, and the system re-evaluated for 


changes in parameters and operating conditions. 


D. Automated State Model Generation Algorithm 


P. C. Krause and Associates (PCK&A) are using a new approach of automatically 
generating state space models in their simulation, and study, of the candidate electrical 
system for the AAAV. This new approach is termed the Automated State Model 
Generation Algorithm [Ref. 4]. As the software is proprietary and not widely available, 
this method will not be used in the simulation and study conducted by this thesis. This 
background is included to ensure a complete and thorough understanding of the methods 


employed by PCK&A in their development and study of the candidate electrical system. 


background is included to ensure a complete and thorough understanding of the methods 
employed by PCK&A in their development and study of the candidate electrical system. 

In the Automated State Model Generation Algorithm the candidate electrical 
system is described by the standard node incidence matrix and elementary branch data 
(e.g. resistances inductances, and back emf's) [Ref. 4]. The standard node incidence 
matrix 1s established by numbering the nodes and branches of the circuit. For example in 
Figure 3-1 below, the node incidence matrix would be of the form shown in Figure 3-2. 
The nodes are labeled vertically on the right while the branch elements are specified 


horizontally across the top. 





Figure 3-1, Example Circuit. [Ref. 4] 
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Figure 3-2, Node Incidence Matrix. [Ref. 4] 


The elementary branch that PCK & A employs to model portions of systems ranging 
from thyristor switches to machine phases is illustrated in Figure 3-3. If the switch is closed 


(state=1), the branch is active. If the switch is open (state=0), the branch is inactive. 


+ У; z 
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ee 
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Figure 3-3, Branch Model. [Ref. 4] 


From this input data the Automated State Model Generation Algorithm will automatically 


establish the composite system state equations, thus relieving the analyst from establishing the 


differential equations by more tedious conventional means [Ref. 4]. Since the state equations are 
now available, linearization and eigenvalue analysis can be performed at the system level. 
Further, the state equations can then be solved numerically using one of a number of well- 


established numerical algorithms (i.e., ones available in ACSL: Runge-Kutta, Gears-Stiff, etc.). 
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IV. ELECTRICAL SYSTEM COMPONENT BACKGROUND 


А. Introduction 


This Chapter includes a description of the selected candidate electrical 
system architecture. In addition it provides background on the major components of the 


selected architecture, specifically the NBC fan drive, ultra capacitor, and batteries. 


B. Candidate Electrical System Description 


Figure 2-6, illustrates a general electrical system architecture for the 
AAAV. From this system architecture PCK&A selected a candidate system, Figure 4-1, 


for initial modeling purposes and to demonstrate the principle of distributed computing. 
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Figure 4-1, AAAV Candidate Electrical System Description. [Ref. 2] 


C. NBC Fan Drive Background 


1. Introduction 


The brushless dc motor is becoming widely used as a small-horsepower control 
motor. The brushless dc motor, however, is really not a dc motor at all; in fact it is a 
permanent-magnet synchronous machine. Nevertheless, when it is supplied from a 
source with frequency equal to the electrical angular speed of its rotor, it becomes a 


brushless dc motor. The nomenclature is adopted not because it looks like a dc motor, but 
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because its torque speed characteristics can be made to resemble those of a dc shunt 
motor with a constant field current. If the permanent-magnet synchronous machine 1s to 
be supplied from a source of frequency that always corresponds to the electrical angular 
speed of its rotor, it is necessary that the rotor position be identified. This can be done in 
several ways, such as with a light beam and a toothed gear, or indirectly by using the 
back emf voltages found by sensing the open-circuit phase voltages, or by using Hall- 
effect sensors which magnetically sense the position of the rotor poles. The stator voltage 
source and rotor electrical angular speed are self-synchronized using a dc-to-ac inverter, 
wherein the transistors are switched on and off based on rotor position. There are 
numerous types of inverters. This thesis will explore the 180? and the 120? voltage source 
inverters. 

In this section a description of the permanent-magnet synchronous machine is set 
forth. Next, the 180° inverter operation is explained. The operation of the 180° inverter 
is then used as a basis for the explanation of the operation of the 120° inverter. This 


explanation is completed by defining a sensorless control method for the 120° inverter. 


2. Description of the Brushless DC Machine 


The term "brushless dc" is somewhat of a misnomer since the brushless de 
machine is actually a permanent-magnet synchronous motor such as the one depicted in 
Figure 4-2. In fact, the only difference between the common notion of a permanent- 
magnet synchronous machine and the brushless dc machine is in the form of the applied 


voltages. 
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In this chapter a background on the brushless dc machine is provided beginning 
with a description of the permanent-magnet synchronous machine in machine variables. 
The machine equations are then transformed into the rotor reference frame. Finally, in 
the context of the rotor reference frame, the form of the applied voltages is specified. 


Additional background is available in references [5 ] and [6]. 


be axis 


as axis 





Figure 4-2, Three-Phase Permanent-Magnet Synchronous Machine. 


a. Machine Equations 


The purpose of this section is to set forth the machine equations for the 
permanent-magnet synchronous machine. The equations presented here are subject to the 


following assumptions: (1) the machine is non-salient, (2) the back emf is sinusoidal, and 
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(3) the machine is wye connected. The stator phase voltage equations describing the 


permanent-magnet synchronous machine are given in vector form by (4.1) 


= 


abcs s’abcs 


= 


y 951 (4.1) 


abcs 


where p is the Heaviside operator d/dt. The flux linkage equations are similarly given by 


(4.2) 


L =M M sin(9, ) 
2 2470 inl 6 а = (4.2) 
-M -M L 


d al 


where A, is a machine parameter dependent on the strength of the permanent magnet. 
Since the machine is wye connected ( the three stator currents must sum to 
zero), (4.2) may be rewritten as (4.3) where L, 1s defined by (4.4). This representation is 


convenient since-the three stator equations become uncoupled. 


sin(@, ) 
Лане =, (ө, - =) (4.3) 


ө, а 
5 
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L.=L+M - (4.4) 


The electromagnetic torque (7?) is established by expressing the co- 


energy, taking the derivative of the co-energy with respect to @,, and multiplying the 


result by the number of pole pairs [Ref. 5]. 
Т, - a c сов(0,)-- E. - i, )sin(, | (4.5) 


b. Machine Equations in the Rotor Reference Frame 


Considerable insight and simplification may be gained by making a 
change of variables which transforms all machine variables (abcs) to rotor reference 


frame variables (qdQs). The transformation from abcs variables to qdOs variables 1s given 


by (4.6) and (4.7) where 7 may represent voltage, current, or flux linkage. 


өле x: I та (4.6) 


cos(@_ ) cos 6 - = cos, + =) 


Re 2 sin(6, ) (ө, à sin, + =) (4.7) 
vea 3 3 
E B E 
p 2 2 | 
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The angle 8, is the rotor electrical angle that necessarily varies at the rotor 


electrical angular velocity. The inverse transformation is 


cos(9, ) sin(6, ) 1 


(к, D = cox 8, - 27. за (4, - 2 | ] . (4.8) 


| cos 6 +=) inl 6 + =) 1 
| 3 3 


Applying the transformation to the machine equations and expanding 


those equations yields 


v; 2 (r, * pL,)i;, * 9, Lir, * 9A, (4.9) 
vi, =(r, + pL, ) iz, - o,L,i;. (4.10) 


зүр\.,., 
m = БЕ (4.1 1) 


It can be seen from (4.11) that if torque is to be constant, then i,, must be 


constant. Hence from (4.9) and (4.10), steady-state operation requires that v, and v;, 


are constant as well. The steady-state line-to-neutral voltages may be given by 
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V. cos(0, + а) 
p T cos, -Z sa] (4.12) 


abcs 
V. cos, + = + “| 


where @ is some arbitrary phase angle. Note, to achieve brushless dc operation, the stator 


voltages must vary at the rotor electrical angular velocity. Then, for any rotor speed, 


(4.13) and (4.14) apply. If the amplitude V, is constant, then so are v;, and vj, . 
vi, =V, cos(a) (4.13) 
v^, — -V, sin(a) (4.14) 


For steady-state conditions, with a —- 0 and o,L, ««r,, (4.13), (4.14), 


(4.9)-(4.11) may be manipulated to yield (4.15) which is identical in form to the steady- 


State torque equation of an armature-controlled dc machine. 


Т, = [ze M (4.15) 
r 


2Д2/к 
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Now that the basic features of brushless dc operation have been set forth, 
the problem of actually obtaining the voltages described by (4.12) is addressed. One 
device used to obtain these voltages in practice is known as an inverter. The inverter is 


the subject of the next section. 


3. Inverter Operation 


There are several methods of achieving stator voltages of the form given by 
(4.12), where the frequency of the voltages is variable. One method is to use a device 
known as an inverter. A typical three-phase inverter is depicted in Figure 4-3. Assuming 
ideal components, the inverter operates without loss and thus a conducting device is 
modeled as a short circuit and a blocking device is modeled as an open circuit. The 
tradeoff in using an inverter is that the voltages produced contain harmonics that increase 
loss in the stator circuits and can introduce torque pulsations. 

There are two types of inverters, the 120° and the 180°. In the 120° case, each 
transistor is on 120° out of a 360° cycle. Likewise, in the 180° case, each oso 15 ОП 
180? out of 360?. Although the 120? inverter-driven brushless dc machine is more 
difficult to handle analytically, it has a practical advantage in that it can be operated 
without any rotor position sensors. In this case, position information is derived indirectly 


from measurement of the open-circuit phase voltages. 
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a. Basic Inverter Operation 


A schematic depiction of the brushless DC inverter system appears in Fig 
4-3. The inverter consists of six transistors, T1-T6, which are used to convert DC power 
to AC power and six free-wheeling diodes, D1-D6, which are necessary to protect the 


transistors from the inductive nature of the load. 





Figure 4-3, Inverter and Machine. 


The first step toward the analysis of the inverter is to apply Kirchhoff's 
Voltage Law (KVL) starting at the negative dc rail (g), proceeding across the lower 
transistor and across the machine phase, and returning to the negative dc rail via the 


neutral point (n). Applying this path for each machine phase yields 
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Varas (4.16) 


у, -УЬ ЖУ, (4.17) 


g 
Veg = Ves + Vag (4.18) 
Kirchhoffs Current Law (KCL) dictates for a wye connection that 
i +i, +i, =0 (4.19) 


Equations (4.19) and (4.3) may be used to show that the sum of the stator flux linkages, 


namely A,,+A,,+4,,,1S zero. Since the sum of the stator currents is zero, and the sum 


cs ? 


of the stator flux linkages is zero, (4.1) shows that 
у БУ РЎ =0 (4.20) 
15 


The neutral to ground (where ground is the lower rail of the dc supply) voltage, v,, , 


obtained by adding (4.16), (4.17), and (4.18), and then subtracting (4.20) which results in 
пе = lv, decr va) (4.21) 
E 
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Substituting (4.21) into (4.16), (4.17), and (4.18) yields the following expressions for the 


phase voltages. 


2 1 1 
Vee = 3 Yag ~ 3 5g 73 « (4.22) 
2 1 1 


= 3 ба T3 Yag Са 4 (4.23) 


2 1 1 
Vie = 27 - 3 Yag - p Vas (4.24) 


b. Six-Step 180? Continuous-Current Inverter 


The logical signals for the switching (commutation) of the transistors are 
shown in Figure 4.4. T; is turned on at 0, = -90° and turned off at O, = +90, at the 
same time T41s turned on. Note, in practice some delay must be present to avoid short- 
circuiting the dc input supply. The gating signals for the top transistors associated with 
the remaining two phases (T3 and Т5) are each phase displaced from the T; signal by 120? 
so that a balanced set of voltages is produced across the load. It is assumed that the turn- 
off time of the transistors is negligible whereupon the transistors become ideal switches. 
At the instant T; is turned off, and T4 1s gated on, the current T, was carrying is diverted 


to the diode D4. The diode continues to conduct until the current decreases to zero. Once 
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the current i,, reverses direction it is carried by T4. This mode of inverter operation is 


generally referred to as the continuous-voltage or continuous-current mode. [Ref. 7] 


мент ккал 


a MAA 
ө Джоан 


0 


5 


Figure 4-4, Transistor Switching Logic for the 180? 
Continuous-Current Inverter. 
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The operation of the 180? inverter 1s set forth in Table 4.1. In the table 
electrical rotor position 1s divided into six 60? segments known as Basic Switching. 
Intervals (BST). The second column specifies the center of each switching interval in 
terms of 0,+¢. The angle @ is known as the advance-firing angle, which allows the 
shifting of the transistor firing sequence relative to the position of the rotor. The third 
column of Table 4.1 defines the transistor firing sequence by indicating which transistors 
are gated on ina given BSI. The transistor firing pattern results in pole (abcg) and phase 
(abcs) voltages which are specified in Table 4.1 [Ref. 7]. Note that once the switching 
pattern is assigned, the phase voltages are completely determined. This fact makes 


programming the 180° inverter rather straightforward in ACSL. 


О,+% Transistors On 


< 
& 





Table 4.1, 180? Inverter Operation. 
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Expressing the phase voltages given in Table 4.1 as a Fourier series in 
electrical rotor position, neglecting all harmonics, and incorporating the transformation to 


the rotor reference frame yields 


pa Ve cos(g) (4.25) 
T 
у” = Ee) (4.26) 
A 


Comparing (4.24) and (4.25) to (4.13) and (4.14), it is apparent that the 180? inverter o 


and @ are equivalent. 


c. 120? Discontinuous-Current Inverter 

The logical signals for the switching (commutation) of the transistors are 
shown in Figure 4-5. T, is turned on at Ө, = 0° and turned off at 9, = +120° ; however, 
T41s not turned on until 60? later at 0, 2180? Hence, the diode in parallel with T4 begins 
conducting when T; is turned off and continues to conduct until the current i,, becomes 


zero. At this point the as winding of the machine becomes open- circuited. The 60?. delay 
in turning T4on is done primarily for the purpose of lessening the possibility of a short 


circuit occurring through T; and T; and to provide an open-circuit period for indirectly 


sensing the rotor position. The voltage v,, in Figure 4-3 is equal to the source voltage 
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v ,, if Ti is turned on or if the diode in parallel with T; is conducting. When T4 1s turned 


on or when the diode in parallel with T4 is conducting, v, is zero. 


Ti | | 
T, | | 
T m m 


16 | | 


0 Vax ж 312 х 2.5 


0 


Y 


Figure 4-5. Transistor Switching Logic for the 120? 
Discontinuous-Current Inverter. 


The operation of a 120? inverter is summarized in Table 4.2, which 15 
much the same as Table 4.1 except for two important differences. The first 1s that only 


two transistors are gated on at any given instant of time. Secondly, question marks 
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appear periodically in the last three columns to indicate that the instantaneous current 


determines the voltages in the corresponding phase. [Ref. 7] 


BSI 


даанн- 


Ө, + 


0° 
60° 
120° 
180° 
240° 
300° 


Transistors On 


1,6 
1,2 
2,3 
3,4 
4,5 
5,6 


Мз 


о Ф о ~ д (9 
а о 


Table 4.2, 120? Inverter Operation. 


The dependence of the phase voltages on the current in the open-circuited 


phase is best illustrated by examining one specific BSI. In BSI II if i,. < 0, then current 


will flow through D3 and thus v,, = V,,. Conversely, if i,,> 0 then current flows | 


through D6 which causes v,, —0. Finally, if i,, — 0, then v,, is the open-circuit voltage 


as derived from (4.1) and (4.3) and given by 


у = 4,0, cos 6, E = 
3 


(4.27) 


If (4.17) is substituted into (4.21), then the neutral to ground voltage is given by 
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Vig = zs ғу, +V cg } (4.28) 


Upon using (4.27) for v,, in (4.28) and then using (4.28) to replace v,, in (4.17), the 
voltage across the lower transistor in leg В (То) 15 


у = Va + 22,0, cos 6, и =) (4.29) 


Similarly, (4.28) can be used to find v,, and v,, from (4.16) and (4.18). 


D L - EE cos, - =) - (4.30) 
2 2 3 
D Б: Vo, и cos 6, E) (4.31) 
2 2 3 


This analysis can be repeated and adapted for each basic switching 


interval (BSI). Table 4.3 reflects the pole voltages, v,,.,, when each phase current, isses» 1S 


abcs ? 


exactly equal to zero. For each phase current, i s, greater than zero, each pole voltage, 


Vaseg » 18 equal to V,.. For each phase current, i 


abcs ? 


less than zero, each pole 
voltage, v,,., , is equal to zero. Table 4.4 similarly reflects the stator phase voltages, 


Vaes >» When the phase currents, i,,.., are equal to zero. 


abcs ? 
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x V, + T cos, + = 
2 2 3 












5 И т 4.0, со5(0, ) 


| У. + 2 24,0, cos, + = 
2 2 3 





b 


Table 4.3, Pole Voltages for BSI I through VI with Non-Conducting 
Phase Open Circuited. 
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a V,- E cos, - | 
2 2 3 






E y, - 1 2,0, cos 6, + = 
2 2 3 





| | 


| 
VI 3009 5,6 
4,0, со5(0,) 


Table 4.4 Stator Phase Voltages for BSI I through VI with 
Non-Conducting Phase Open Circuited. 


ai V,- B cos 6, E 
2 2 3 
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d. Sensorless Control 


Brushless dc motors are usually controlled using dc-ac inverters as 
previously discussed. The motor requires a rotor position sensor for starting and for 
providing the proper commutation sequence to turn on the power devices in the inverter. 
Hall effect sensors can be used; however, they must be specially mounted in the interior 
of the stator. Further these sensors are temperature sensitive, limiting the operation of the 
motor to below about 75°C. Due to these limitations of the motor operation with Hall 
devices, sensorless operation of the PM brushless motor is receiving more attention. 

Several schemes for position sensorless operation of PM motors have been 
reported and reviewed. They are based on the following techniques. [Ref. 8] 

l. Position sensing using the back emf of the motor. 

2. Position information using the stator third harmonic component. 

3. Position information based on detection of the conducting interval of 
free-wheeling diodes connected in anti-parallel with the power 
devices. 

4. Position information using the motor parameters, and the measured 
terminal voltages and currents. 

Position sensing using the back emf of the motor is being utilized for 

control of the NBC fan motor. The sensed back emf can be used in several ways to 


obtain the proper switching sequence for the inverter power devices. Some or these 


methods are: 
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e The zero-crossing approach based on detecting the instant at which the 
back emf in the unexcited phase crosses zero or reaches a predetermined 
level. 

e Phase-lock-loop (PLL) technique to lock on to the back emf waveform in 
the unexcited phase winding during each 60 degree excitation interval. 

e Back emf integration approach, which provides a commutation pulse when 
the integrated rectified back emf waveform reaches a preset threshold 
value. 

The zero-crossing approach method is utilized in the NBC fan drive and therefore is 
discussed. In a 120? inverter-driven permanent-magnet brushless dc machine only m of 
the three phases are excited at any one time. The back emf voltage in the unexcited phase 
can be measured to establish a switching sequence for commutation of the power devices 
in the inverter. Based on the rotor position, the power devices are commutated 
sequentially every 60 degrees to continually synchronize the phase excitation with the 


rotating rotor magnetic field. 

The pole voltages v,, , v,, , v, , and the neutral voltage v,, are measured. 
At the instants of the zero crossing of the back emf waveform, the pole voltage is equal to 
Va 


С 


. In 





the neutral voltage. When the back emf of phase A crosses zero, then v,, =v,, = 5 


order to use this zero-crossing point to derive the switching sequence, this point has to be 


phase shifted by 90 degrees [Ref. 8]. 
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D. Ultra Capacitor Background 


As shown in the candidate electrical system description, an ultra capacitor bank 
will be modeled and validated. This ultra capacitor may be incorporated in the AAAV 
electrical system in an effort to capture regenerative energy. Regenerative energy could 
be produced as a result of a sudden impact on the vehicle's main gun turret. This energy 
could be forced back into the high power distribution bus thus potentially ет the 
system. Theoretically, the ultra capacitor bank will absorb this energy and prevent any 
damage. 

Ultra capacitors are being considered over conventional capacitors due to their 
high capacitance relative to their small compact size. Since the technology 15 new, a short 
discussion on ultra capacitors will follow. 

Conventional capacitors store energy in the form of separated electric charge. The 
greater the area for storing charge and the closer the separated charges results in greater 
capacitance. A conventional capacitor achieves its cross-sectional area from plates of flat, 
conductive material. To realize high capacitance, this material can be wound in great 
lengths to increase its surface area. A conventional capacitor separates its charged plates 
with a dielectric material, sometimes plastic or paper film, or even a ceramic. These 
dielectrics can be made only as thin as the available films or other applied material. 

An ultra capacitor gets its surface area from a porous carbon-based electrode 
material. The porous structure of this material allows its surface area to approach 2000 
square meters per gram, much greater than can be accomplished using flat or textured 


films and plates. An ultra capacitor’s charge separation distance is determined by the size 
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of the ions in the electrolyte, which are attracted to the charged electrode. This charge 
separation is much smaller than can be accomplished using conventional dielectric - 
materials. 

The combination of enormous surface area and extremely small charge separation 
gives the ultra capacitor its outstanding capacitance (1 Farad, 100 Volts) relative to 


conventional capacitors. 


Е Battery Background 


As shown in the candidate electrical system description, two battery banks, one 
connected to the Vetronics bus and the other to the auxiliary bus (see Figure 4-1), will be 
modeled and validated. Therefore, a brief discussion on basic battery concepts is 
included. 

A battery is a device that converts the chemical energy contained in its active 
materials directly into electric energy by means of an electrochemical oxidation-reduction 
(redox) reaction. This type of reaction involves the transfer of electrons from one material 
to another through an electric circuit. In a non-electrochemical redox reaction, such as 
rusting or burning, the transfer of electrons occurs directly and only heat is involved. 

While the term battery is often used, the basic electrochemical unit being referred 
to 1s the cell. A battery consists of one or more cells, connected in series or parallel, or 
both, depending on the desired output voltage or capacity. 


The cell consists of three major components. 
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1. The anode or negative electrode — which gives up electrons to the external 
circuit and is oxidized during the electrochemical reaction. 

2. The cathode or positive electrode — which accepts electrons from the external 
circuit and is reduced during the electrochemical reaction. 

3. The electrolyte — the ionic conductor — which provides the medium for 
transfer of electrons, as ions, inside the cell between the anode and the 
cathode. The electrolyte is typically a liquid, such as water or other solvent, 
with dissolved salts, acids, or alkalis to impart ionic conductivity. 

The lead-acid battery, as will be used in the AAAV, uses lead dioxide as the 

active material for the positive electrode and metallic lead, in a high surface area porous 
structure, as the negative active material. The electrolyte is a sulfuric acid solution 


[Ref. 9]. 
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M. COMPONENT MODELING AND VALIDATION 


A: Introduction 


This chapter contains a discussion of the modeling efforts applied to the NBC fan 
drive, ultra capacitor, and battery unit. It also provides results and conclusions for each 


model. 


B. NBC Fan Drive Model 


The circuit diagram of the NBC fan motor and its associated three-phase 
inverter drive is shown in Figure 5-1. The resistance and the inductance of the cable 
connecting the motor drive to the high power 28-volt bus are included in the diagram as 


Yline and Liine. The input capacitance of the inverter is given by Ca. 


2“ 
—> 


+ Thins Lis 


Vip Ун be Ce a 
T, « 


Figure 5-1, NBC Motor Drive System Circuit Diagram. 
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As discussed in Chapter 4, the equations describing the brushless dc motor are expressed 


as. 
уг, = (г, + pL,)ir * o, Ljis * 0,4; (5.1) 
vi, =(r, * pL,) is —0,L, iz, (5.2) 
T, = poz (5.3) 
where 
@„= о (5.4) 


апа Р is the number of poles. These equations together with the transformation of the 
phase voltages into the rotor reference frame and the inverse transformation of the rotor 
reference frame currents back to machine variables are coded in ACSL in a 
straightforward manner. 


The mechanical speed is established by integrating 


po, 7 (T, - T, J, | (5.5) 
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where T, is the electromagnetic torque, T; is the load torque, and Jz is the inertia of the 
motor and the connected fan load. It is assumed that the load torque is proportional to the 


square of the rotor speed. In particular, 


na Kk o (5.6) 


m 


The actual inverter is a 120? discontinuous-current inverter based on sensorless 
control. The current model being exercised by PCK&A is that of the NBC fan motor 
with a 180? continuous-current inverter using Hall-effect sensors to detect rotor position. 
The first step in developing an accurate model was to build a model similar to that being 
currently exercised. This model would allow for a better understanding of the inverter 
and motor assembly. 

The 180? inverter model was developed and exercised. The transistor gating 
signals are derived based on the waveforms illustrated in Figure 4-4. The phase voltages 
are then determined from (4.22) — (4.24). The current drawn by the inverter is 
established by summing the phase currents corresponding to the gated transistors 
connected to the top rail. The parameters used are shown in Table 5.1. The parameters 
were derived using manufacturer provided data. The complete ACSL model can be 
found in Appendix A. For convenience, a fixed 28-volt ideal source was used as the 


input to the inverter. 
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1, 0.005 mH Stator self Calculated from manufacturer 
І Inductance supplied data 
BOUT mmm (average of upper and lower values) 
r 0.003 Q Stator resistance | Calculated from manufacturer 
| Supplied data 
P | I d (average of upper and lower values) 
0.00344 V-s/rad Back emf Calculated from manufacturer. 
E ses LU eed IL 


Number of poles | Communication from manufacturer 


y 1.171Х10” Котп Inertia Estimated 


1.128X10~ Load constant Estimated from manufacturer 
Supplied data 


Table 5.1, NBC Fan Motor Drive Parameters. 



















Figure 5-2 shows the first few cycles of the stator currents, jas, ips and ics, Upon 
Startup of the 180° inverter model using the data derived from communication with the 
manufacturer. With the low value of stator resistance, these currents illustrate a large 
startup transient, on the order of 4000A for ia. Since rated current for the drive is 
approximately 125A, this current draw represents a significant transient load on the high- 


power bus. 
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Figure 5-2, Initial Stator Currents of NBC Fan Model 
with 180? Continuous-Current Inverter Upon Startup. 


Figure 5-3 shows simulation results for the drive extending from startup to steady 
state. The variables plotted include @, — rotor speed in rads/s, 7; — electromagnetic 


torque in Newton-meters, i4, - stator current in Amps, and ij; — inverter input current in 


Amps. The large spike in inverter input current is readily apparent as are the considerable 


S 


current harmonics resulting from the applied six-step voltage waveforms. The motor 


speed settles out at a final value of about 1250 rad/sec. 
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Figure 5-3, Plot of Start-up of NBC Fan Model with 
180? Continuous-Current Inverter. 
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Figure 5-4 illustrates steady-state waveforms for the 180° inverter drive. In 
particular, the figure shows variables v,, - stator voltage in Volts, i4 - stator current in 
Amps, and ij; — inverter current in Amps, 7; — electromagnetic torque in Newton-meters, 
and o,, — rotor speed in rads/s. The 6c, dominant harmonic in ij; is apparent and the 


peak value of i4 appears to be on the order of 800A. This 800A amplitude is much 
larger than the anticipated 125A, though the final speed of the device is rather close to the 


anticipated rated speed (1310 rad/sec). The difference may be due to the drive actually 


being a 120° inverter drive or to an inaccurate assumption on the motor parameters. 
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Figure 5-4, Plot of Steady State of NBC Fan Model with 


180? Continuous-Current Inverter. 


Research then began on implementation of a 120? discontinuous-current inverter 
using a sensorless control scheme. As discussed in Chapter IV, paragraph C, subsection 
3.b, sensorless control can be initiated by indirectly assessing the rotor position using the 
back emf of the motor. The back emf can be monitored to detect the instant at which the 
back emf in the unexcited phase crosses zero or reaches a predetermined level. The 
instant at which this happens is then used to establish a switching sequence for 
commutation of the transistors in the inverter. 

The main nuance between the simulation of the 120° and the 180° inverter is 
establishing which devices are conducting. In particular, when neither transistor in the 
inverter leg is gated, then one of the diodes may be conducting or the phase could be 
open-circuited. In addition depending on the back emf, the open-circuited phase voltage 
may change and cause one of the diodes to turn on. The logic used in the simulation for 


establishing the status of D; and D, is illustrated below. 


IF ((T1) .OR. (T4)) THEN 


D1 = false. 
D4 = „false. 


ELSEIF (ias .GT. hold) THEN 


D4 = .true. 
D1 = . false. 


ELSEIF (ibs .LT. -hold) THEN 


D1 = .true. 
D4 = „false. 


ELSEIF ((Vagoc .LT. (Vdc + eps)) .OR. (Vagoc .GT. -eps)) THEN 


ПІ = .false. 
D4 = .false. 


56 


ELSEIF (Vagoc .GT. (Vdc + eps)) THEN 
Di mue. 

ELSE 
D4 = .true. 


ENDIF 


Note, if neither transistor is gated and the phase current is positive, then D4 must conduct. 
Conversely, if the current is negative, the Dı must conduct. Otherwise, the back emf 
must be assessed. If it is greater than the input voltage (Vac), then the upper diode is 
turned on. If it is less than zero volts, then D4 is turned on. If the open-circuit voltage 
remains between zero volts and Vgc, then the phase remains open circuited. 

The complete ACSL code for the NBC Fan with the 120° discontinuous-current 
inverter can be found in Appendix B. The 120° inverter model was developed and 
exercised. Тһе A4; parameter was adjusted through simulation to yield steady-state stator 
current waveforms consistent with manufacturer data. The complete set of parameters is 


documented in Table 5.2. Once again a fixed 28-volt source was used as the input to the 


inverter. 


w 


I 0.005 mH Stator self Calculated from manufacturer 
у Inductance supplied data 
(average of upper and lower values) 


ғ, 0.003 О Statorresistance | Calculated from manufacturer 
Supplied data 
(average of upper and lower values) 


0.006 V-s/rad Back emf Adjusted iiie. — E simulation 
constant 


Number of poles | Communication — 73 manufacturer 


J, 1.171Х10” Кетп Inertia Estimated 


1.128X10” Load constant Estimated from manufacturer 
supplied data 


Table 5.2, NBC Fan Motor Drive Parameters. 































Figure 5-5 shows the first few cycles of the stator currents, ias, ibs, and ies, upon 
startup of the 120? inverter model assuming that rotor position is available. Similar to 


the 180° model, the currents illustrate a large starting transient, reaching levels of | 


approximately 2600A. Also, following the initial swing, the currents assume the classical 


120* waveform, where the open-circuit waveform is apparent. 
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Figure 5-5, NBC Fan Model Plot with 120? 
Discontinuous-Current Inverter. 


Figure 5-6 illustrates a complete startup to steady state and documents 
Q,,- rotor speed in rads/s, T, — electromagnetic torque in Newton-meters, ia, - stator 


current in Amps, and ijj — inverter current in Amps. 
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Figure 5-6, Plot of Start-up of NBC Fan Model with 
120? Discontinuous-Current Inverter. 
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Figure 5-7 shows two cycles of the steady-state variables: va - stator voltage in 
Volts, i4 - stator current in Amps, and ij; — inverter current in Amps, 7, — 


electromagnetic torque in Newton-meters, and @,, — rotor speed in rads/s. The inverter 


current once again exhibits a6q@, fundamental variation, but it clearly has a different 


harmonic content than in the 180° model. The a-phase terminal voltage closely 
resembles the back emf as it should with the quasi-rectangular phase current waveforms. 
The a-phase stator current has a peak value of about 120A, inline with what the 


manufacturer data predicted. 
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Figure 5-7, Plot of the Steady State of NBC Fan Model with 
120? Discontinuous-Current Inverter. 


The simulation of the NBC fan model with 180? continuous-current inverter was 


run again using the parameters in Table 5.2 to illustrate the effect of varying A7. The 
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results can be found in Appendix C. Also, the simulation of the NBC fan model with 
120? discontinuous-current inverter was run again using the parameters in Table 5.1 to 
illustrate the effect of а A}. The results can be found in Appendix D. 

After successful implementation of the NBC fan model with 120° discontinuous- 
current inverter with a known rotor position, work began on implementing sensorless 
control of the motor. In order to accomplish this the open-circuit voltage signal would 
need to be processed into a signal that could be utilized to gate transistors T; through Tg 
as discussed in Chapter IV, paragraph C, subsection 3d. Figure 5-8 depicts the signal 
processing block diagram that was implemented in order to process the open-circuit 


voltage. 


ае Хас Xy Усо 
— 7 G6/) б) |> б) ү ————>э 


Figure 5-8, Open-Circuit Voltage Signal Processing Block Diagram. 


0; 





а (5) = 2! (5-7) 
6, (5)= = (5-8) 
5 
TS 
G,(s)= aa (5-9) 
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where œw, = 10,000 , a = 2500, r= 0.004 . 

G;, a low-pass filter, smoothes out the choppy terminal voltage waveform (vais). 
G5, an integrator, introduces a 90° phase shift. G; serves to block any de component that 
may be introduced when the integrator starts. The zero-crossings of veo (together with 
analogous Vg, and vs signals from phase A and B) may be used for gating 7;— 7¢. Figure 
5-9 documents the start-up to steady-state waveforms from the block diagram with the 
machine operating with rotor position feedback. In Figure 5-9 the integrator is disabled 
until near the first peak of vgig- . This ensures that x, does not have a large dc offset that 
would then need to settle out. Figure 5-10 documents the steady-state waveforms from 
the block diagram with the machine operating with rotor position feedback. The 90° 


phase shift between vai (and xj.) and xpe (and veo) is readily apparent. 
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Figure 5-9, Plot of the Start-Up of NBC Fan Model with 
120? Discontinuous-Current Inverter. 
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The logical statements required to transition to sensorless operation are illustrated 


below. 


IF ((vco .LT. 0.0) .AND. (vbo .LT. 0.0)) THEN 


Tla=.FALSE. 
T2a = . TRUE. 
T3a = . TRUE. 
T4a 7» .FALSE. 
T5a = .FALSE. 
T6a = .FALSE. 


ELSEIF ((vbo .GT. 0.0) .AND. (vao .GT. 0.0)) THEN 


Tla =.FALSE. 
Т2а = .FALSE. 
ТЗа = .TRUE. 
T4a = .TRUE. 
T5a = .FALSE. 
T6a = .FALSE 


ELSEIF ((vao .LT. 0.0) .AND. (vco .LT. 0.0)) THEN 


Tla = FALSE. 


T24 SE 


T3a = .FALSE. 
T4a = .TRUE. 
T5a = .TRUE. 


T6a = .FALSE. 


ELSEIF ((vco .GT. 0.0) 


Tla = .FALSE. 
T2a = .FALSE. 
T3a = .FALSE. 
T4a = .FALSE. 


T5a = . TRUE. 
Тба = ТЕСЕ. 


ELSEIF ((vbo .LT. 0.0) 


Tla=.TRUE. 


T2a = .FALSE. 
T3a = .FALSE. 
T4a = .FALSE. 
T5a = .FALSE. 


Тба = .TRUE. 


ELSEIF ((vao .GT. 0.0) 


ENDIF 


Tla=.TRUE. 
T2a = .TRUE. 


T3a = .FALSE. 
T4a = .FALSE. 
ба =- FALSE 
T6a = .FALSE. 


.AND. (vbo .GT. 0.0)) THEN 


AND. (vao .LT. 0.0)) THEN 


.AND. (vco .GT. 0.0)) THEN 


Details for the entire simulation are documented in Appendix B. 
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Figure 5-10, Steady-State Waveforms From Block Diagram with the 
Machine Operating with Rotor Position Feedback. 
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The following simulation study was conducted to verify the feasibility of 
sensorless operation. Rotor position information was used to start the machine (no back- 
emf voltages are initially available). At t=0.05, Figures 5-11 through 5-12 illustrate the 
signals derived from the actual rotor position information (Tj, through T¢,). Clearly the 
signals T,, through T¢, are very nearly in phase with T, through T¢, respectively, and can 
be used to gate the transistors. Figure 5-14 depicts transitioning from using the he 
feedback signals to using the sensorless derived signals. The rotor maintains 


synchronism, with only a small drop in speed. 
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Figure 5-11, Plot of the Switching Signals of NBC Fan Model with 
120? Discontinuous-Current Inverter. 
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Figure 5-12, Plot of the Switching Signals of NBC Fan Model with 
120? Discontinuous-Current Inverter. 
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Figure 5-13, Plot of the Switching Signals of NBC Fan Model with 
120° Discontinuous-Current Inverter. 
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Figure 5-14, Plot of the Rotor Feedback Signals and 
the Sensorless Derived Signals. 


While the model introduced does implement sensorless control with good results, 
the problem of transitioning to sensorless control from rotor position feedback must still 


be resolved. This transition is an area that can be addressed in future wok. 
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С: Ultra Capacitor Model 
Pinnacle Research Institute, Inc. provided tabular data on their 100V, 1F, 5000J 
Ultra-Capacitor. The manufacturer's measured data was entered into a MATLAB file, 
ultracap.m. This file is included as Appendix C. The following plot illustrates the 


manufacturer's measured impedance versus frequency data. 


Impedance vs. Frequency 





15. | 


Impedance (Ohms) 





10? 10! 107 10° 
Frequency (Hz) 


Figure 5-15, Ultra Capacitor Impedance vs. Frequency 
(Manufacturer’s Measured Data). 


A computer model of the ultra capacitor would need to produce comparable 
results when exercised over the same frequency range. An attempt to duplicate the actual 


results would be initiated using an equivalent series resistance (ESR) as in Figure 5-16. 
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| SS 
Figure 5-16, ESR Circuit Diagram. 


The input impedance of the circuit is found by inspection to be 


1 
Z =| R+— 5.10 
eq | = ( ) 
which can be rewritten as 
4з + zc 
De ps AE (5:1 


An initial guess for the constants R and 1/RC was formulated using the PRI, Inc. 
frequency response data. After many iterations of manipulating these constants the 


following plot of impedance vs. frequency was achieved. 
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Impedance vs. Frequency 
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Figure 5-17, Ultra Capacitor Impedance vs. Frequency (Model). 


Figure 5-18 illustrates both actual data points and model generated data and 
confirms that the model is valid as a reasonable representation of the actual ultra 


capacitor. 
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Figure 5-18, Ultra Capacitor Impedance vs. Frequency (Actual and Model). 


The constants R and D were found to be: 


R- 0.045 ohms - (5.12) 
1 E 
——- 535s (5.13) 
RC 


and the transfer function representation for the ultra capacitor would be characterized as 


- 


_ 0.045(S + 5.35) 


iy < 


(5.14) 
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D. Battery Model 


The AAAV electrical Jide consists of two battery banks, one connected to the 
Vetronics bus and the other to the Auxiliary bus. The bank connected to the Vetronics 
bus consists of two series connected 12-Volt (nominal) batteries. The bank connected to 
the auxiliary bus consists of three parallel branches with two series-connected 12-Volt 
batteries per branch. In the current Purdue Distributed Computing AAAV system 
representation, the batteries are modeled as a Thevenin equivalent circuit, where the 
open-circuit voltage and the Thevenin resistance can be adjusted to correspond to the 
chemical and thermodynamic state of the battery. The resulting circuit diagram is shown 
in Figure 5-19. The associated parameters are summarized in Table 5.3. In Figure 5-19, 
51, 52, апа 5з аге closed for starting. After start-up S? and S3 are opened, and S, is left 
closed during operation. 

When no battery 1s connected to the Vetronics bus, a bus voltage oscillation was 
observed. The stiffening action of the battery seems to mitigate this instability. It is 
hypothesized that the oscillation is due to a limit cycle imposed by the generator voltage 


regulator, particularly at low power. 
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Figure 5-19, Circuit Diagram of Battery Banks. 












Ro; Vetronics bus battery series | 0.007648 О Calculated using 
ios ЕКЕ 
Ri Auxiliary bus battery series | 0.00255 OQ Calculated using 
ENP ANN PNE 
Е Resistance of cable between | 0.0002 Q Estimated 
A O _ 
A memi 
inductance 
MEM arch Na 
inductance 


Table 5.3, Battery Parameters. 
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The model of Figure 5.19 is preliminary and further laboratory testing has been 
conducted on the 800U Optima battery. Of particular interest was the battery's frequency 
response at various levels of charge. The circuit shown in Figure 5.20 was constructed 
and the parameters for the test circuit are summarized in Table 5.4. The purpose of К 15 
to establish a quiescent load condition on the battery. Transistor Qj is gated at the user- 
specified frequency, with a 50% duty cycle. With R; larger than Rı, switching Qı results 
іп a smaller rectangular sink of current from the battery. Therefore, item consists of a 
steady-state dc level with a small rectangular ripple. This AC component of current flows 
through the battery resulting in an AC component of voltage across the battery (Vterm). 
Therefore, if the assumed internal battery model is that given by Figure 5-21, the AC 
component of the measured voltage will facilitate the calculation of the impedance given 
by Iban and Las. The current was monitored with a Tektronix, 8-Amp, clamp-on, current 
probe, while the terminal voltage was monitored with a direct connection to the 
oscilloscope and with a digital handheld multi-meter. Both signals were displayed on a 
Tektronix, 60 MHz, digital storage oscilloscope. The ratio of the amplitudes of the 
ripples in Vterm and item provides a measure of the impedance magnitude, while the phase 
angle (0) between those signals relates how the impedance divides between real and 


reactive components. Figure 5-21 reflects the assumed internal battery representation. 
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Figure 5.20, Circuit Diagram for Optima Battery Frequency Response Test. 


Th att Ly att 1 term 
nn _ > 
F 
Yes 


Figure 5-21, Assumed Internal Battery Representation. 






Symbol | Description Measured | General Information 
Value 
КЖ Resistor 4 CEN O 4X 16Q resistors in parallel, 
10 W each 


= din dii Resistor im 10 (1) 39 €) 


E | Power MOSFET (2N6796) к 8.0 A, 100 V 
к Optima 800U Battery 20 56 A-hr 


Table 5.4, Test Circuit Parameters. 










Testing was conducted at frequencies of 20 Hz to 20 kHz. The upper end of 
20kHz was deemed sufficient since the ripple frequency out of the 12-pulse 
generator/rectifier is less than 20 kHz. The testing was conducted at three different charge 
states, (100%, 75%, and 50%). The charge state of the battery is determined using the 
manufacturer supplied capacity of 56 A-hr. Full battery specifications are provided in the 


Manufacturer Data Sheet listed in Appendix F. Given that 


1 A—hr «15 X 3600 5 = 3600 С, (5.15) 
5 | 


the stored charge (qstoreg) iS determined as 





9 она = 29 A-—hr X 3600 = 201,600 С . (5.16) 


hr 
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Since charge is related to current by 


4- [iat | (5-17) 


then for constant current (7), equation (5-17) becomes 


а=1 М (5.18) 


By drawing a fixed current (/) over a specific time interval (At), the amount of charge left 
in the battery can be estimated. 
Initial test results of the battery at full charge, V = 12.8 Volts, I =3.5A, and f= 


2kHz are depicted in Figure 5-22. 
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Figure 5-22, Test Results with Battery at Full Charge, 
V = 12.8 Volts, I =3.5A, and f = 2kHz 
СНІ = current (1 Amp/div), CH2 = voltage (1 Volt/div). 


The plot reveals V (mVpp) = 3.0, I (mApp) = 310,9 (тай), = п, giving an internal 
resistance, R (mQ) = 9.67, and internal inductance, L (mH) = 0.0. 

The battery was then drained at a current level of approximately 15.5A for 52 
minutes to arrive at a charge state of about 75% rated. At 75% charge test results of the 


battery with V = 11.9 Volts, I =3.3A, and J = 2kHz are depicted in Figure 5-23. 
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Figure 5-23, Test Results with Battery at 7596 Charge, 
V — 11.9 Volts, I 23.3A, and f » 2kHz 
CH1 = current (1 Amp/div), CH2 = voltage (1 Volt/div). 


The plot reveals V (mV pp) = 2.5; I (mA pp) = 300, © (rads), = 7, giving an internal 


resistance, R (mQ) = 8.33, and internal inductance, L (mH) = 0.0. 


The battery was then drained at a current level of approximately 15A for 55 
minutes to arrive at a charge state of about 50% rated. At 50% charge test results of the 


battery with, V = 11.8 Volts, I =3.25A, and f= 2kHz are depicted in Figure 5-24. 
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Figure 5-24, Test Results with Battery at 50% Charge, 
V = 11.8 Volts, 1 =3.25A, and f = 2kHz 
CHI = current (1Amp/div), CH2 = voltage ( 1 Volt/div). 


The plot reveals V (mV pp) = 2.7, I (mApp) = 290, O (rads), = x, giving an internal 
resistance, R (mQ) = 9.31, and internal inductance, L (mH) = 0.0. 


Complete laboratory test results for all frequencies are summarized in Table 5.5. 
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100% Charge 75% Charge 50% Charge 


a 





Table 5.5, Laboratory Test Results. 


Test results conclude that the internal resistance of the battery does exhibit a small 
change as a result of frequency and charge-state of the battery. Further, tests conclude 
that the internal inductance of the battery is negligibly small. It appears reasonable that 


the battery can be adequately represented by an equation of the form 


Vierm = Vint ~ ^ bau Dos (БЕШЕ) 
If Vem 1S the same as the bus voltage and it is a state variable then 
Vua — V : 
ee — int bus ( 5-20) 
Frau 


where i,,,, 1s then available to evaluate the derivative of the bus voltage. 
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VI. CONCLUSION 


A. Synopsis of Work Completed 


las Introduction 


The Advanced Amphibious Assault Vehicle (AAAV) was introduced 
and a general system description was set forth. The introduction and description are 
critical to understanding where and how the AAAV fits into the Marine Corps 
amphibious triad (AAAV, V-22, LCAC) to execute the concepts of Operational 
Maneuver from the Sea and Ship to Objective Maneuver. The AAAV”s operating modes 


and associated mode characteristics were also documented and developed. 


2 Electrical System Architecture 


The entire electrical system architecture is introduced in block diagram 
form and the system's primary components are listed and discussed. A discussion on the 
theory behind distributed computing techniques through networking is included. This 
technique is the foundation of P. C. Krause and Associates’ troubleshooting tool that can 
be used by the Marine Corps in addressing current electrical system issues that must be 
resolved. An Advanced Continuous Simulation Language (ACSL) background is : 
included. This software is used in this research to model and simulate various 


components of the electrical system. 
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S Candidate Electrical System Architecture 


From the entire electrical system architecture, a scaled down candidate 
system architecture 1s developed. Components from the candidate system are tested, 
modeled, simulated, and validated. The Nuclear, Biological, and Chemical (NBC) 
fan/filter motor drive, an ultra capacitor, and a battery unit are the components selected 


for research in this document. 


4. NBC Fan/Filter Motor Drive 


A simulation of the NBC fan motor drive with a 180° continuous- 
current inverter using Hall-effect sensors to detect rotor position was modeled in ACSL. 
The simulation allowed for a thorough understanding of the inverter and motor assembly, 
and also served as a tool for gaining an understanding of ACSL simulations and 
representative signal waveforms. 

Efforts then shifted to implementing a NBC fan motor drive model with the 120° 
discontinuous-current inverter based on sensorless control. This model was achieved in 
two stages. First, the 120° discontinuous-current inverter model was developed 
assuming the availability of rotor position. The second step was to incorporate sensorless 
control into the simulation. This was achieved by indirectly sensing the rotor position 
using the back emf of the motor. The switching sequence for commutation of the 
transistors in the inverter was established using this back emf. The results were 
consistent with theoretical expectations ensuring this model was a valid representation of 


the NBC fan motor and inverter assembly. 
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5. Ultra Capacitor 


Tabular data was collected from Pinnacle Research Institute Inc., on 
the 100V, 1F, 5000J, ultra capacitor. This data was analyzed and plotted using MATLAB 
and an impedance vs. frequency plot was developed. An attempt to duplicate this plot 
using a model with an equivalent series resistance (ESR) was initiated. A transfer 
function was developed and an initial guess was made for constants. By manipulating the 
transfer function constants, an impedance vs. frequency plot matching that plot from the 


manufacturer's data was achieved. 


6. Battery Unit 


The area of concern for the battery unit was its frequency response 
characteristics at different charge states. A test circuit was constructed that would load 
the battery drawing a constant DC level of current and an AC component of current. The 
current and the voltage levels were displayed on an oscilloscope and monitored. From 
the steady-state and ripple current and voltage waveforms, an internal impedance was 
calculated. These tests were conducted at three charge states (100%, 75%, and 50%) of 
the battery over a frequency range of 20 Hz to 20 kHz. The upper frequency level of 20 
kHz was based on the ripple frequency of the 12-pulse generator/rectifier. Test results 
concluded that the internal impedance of the battery was almost purely resistive with a 
negligibly small internal inductance. Further, tests concluded that the internal resistance 


of the battery did fluctuate as a result of charge state and frequency. 
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B. DOD Relevance 


The AAAV will join the Marine Corps inventory to complete the amphibious 
triad which also includes TM V-22 Osprey, tilt-rotor aircraft, and the Landing Craft, Air 
Cushion (LCAC). The AAAV will provide the principal means of water mobility, land 
mobility, and direct fire support to Marine Corps infantry units. 

As the AAAV travels through its prototyping stage in an effort to successfully 
pass its operational testing, electrical system issues and problems have arisen which must 
be addressed before entering into full-scale production. Computer simulation of the 
electrical system will facilitate the analysis and evaluation of hardware design prior to 
and in concert with actual production. Accurate simulation will provide invaluable 
insight into design flaws or system limitations before extensive capital investments are 


made on hardware. 


С. Further Work 


Complete simulation for the AAAV’s entire electrical system is the next step in 
providing the Marine Corps with a tool where by it can successfully analyze and СІНЕ 
the system. Once a complete and accurate representative model exists, a distributed 
computing network must be established in order to exercise the entire model. Funding 
for such a network at NPS has not yet been obtained and a future researcher may be able 
to pursue possible avenues for funding. 


Other research that can begin now includes a paralleling strategy for the auxiliary 


generator. The current auxiliary generator is incapable of supporting predicted power 
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requirements. Finally, research could be initiated in an attempt to resolve instability 
problems in the electrical system during low power operation. General Dynamic's 


engineers recently discovered low power instability problems during prototype testing. 


9] 
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APPENDIX A. [ACSL CODE FOR NBC FAN MODEL WITH 180" INVERTER] 


PROGRAM NBC Fan Model w/ 180 Degree Inverter 


INITIAL 


MAXTERVAL maxt = 1.0e-6 
MINTERVAL mint = 1.0e-7 
CINTERVAL cint = 1.0е-4 
ALGORITHM ialg = 5 


CONSTANT tstop = 0.165 


"_--initial circuit constants---" 


CONSTANT rL=0.01 
CONSTANT LL=.000001 
CONSTANT Cdc=.036 

' CONSTANT Vdc-28.0 


LOGICAL SA 
LOGICAL SB 

LOGICAL SC 

"----switch initial conditions----" 
SA = .false. 

58 = .2а15е. 

SC 3 .£alse. 

RND uer аз ы ны 


DYNAMIC 


TERMT (T GZ. tston) 


DERIVATIVE 
l"Swiees States: Gene =ate or” 


CONSTANT 92 
sopa = 2.0 


Ul 


21412925 
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PROCEDURAL (thetar = thetarm) 


thetar = (Poles/2.) * thetarm 
(thetar-twopi) 


(thetar) 


SCHEDULE reset .XP. 
SCHEDULE reset2 .XN. 


END ! End of Procedural 


I^Firing of Switches" 


IF ((thetar .GE. 0.0) .AND. (thetar LT. 5057) 
SA - .true. 
SB = .false. 
SC = .false. 
ELSEIF ((thetar .GE. pi/6.) .AND. 
SA - .true. 
SB = .true. 
SC Басе; 
ELSEIF ((thetar .GE. (pi/2.)) .AND. (thetar .LT. 
SA = .false. 
SB = .true. 
sc = .false. 
ELSEIF ((thetar .GE. (5.*pi/6.)) .AND. (thetar 
SA = .false. 
SB = .true. 
SC =E true. 
ELSEIF ((thetar .GE. (7.*pi/6.)) .AND. (thetar 
SA = .false. 
SB = .false. 
SC = .true. 
ELSEIF ((thetar .GE. (3.*p1/2.)) .AND. (thetar 
SA - .crue. 
SB = .false. 
SE =-.ETUe. 
ELSEIF ((Chetaz GE. (ІТ 79170000 замра (пева 
SA = .true. 
SB = .false. 
$С = га зе: 
ENDIF 


!"Assume lower switches are conducting" 
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THEN 


(thetar Er 2972/22 THEN 


(SA Da THEN 
ГОТ PL Sm I E THENG 
.LT. (3.*pi/2.))) THEN 
.LT. (11.*pi/6.))) THEN 

CUT. ME THEN 


PROCEDURAL (Vas, Vbs, Ves, 111 = SA SD. SC, Vdoc*? "3s 2215s. ісе 
Vap = 0.0 
iia = 0.0 
Vbp = 0.0 
iib = 0.0 
Vep 059 
Lie Saor 


!* Check if top switch is closed" 


IF (SA) THEN 


Мар = Мае 
lia = las 
ENDIF 
IF (SB) THEN 
Vbp = Уас 
iib = ibs 
ENDIF 
IE (SC) THEN 
Vcp = Vdc 
iic = ics 


ENDIF 


!" Establish the phase voltages" 


Мар = (Vap + Vbp + Vep)/3.0 
Vas = Vap - Vnp 
Vbs = Vbp - Vnp 
Ves = Vep - Упр 


t" Find the inverter input current" 
Liz = ila + 11b + lie 
END ! End of Procedural 


Pe Link Dynamic Equations" 


weap = 28.0 


T Brushless DC Motor Constants” 
Ls - .000005 
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“es = 0.003 
lamref = 0.006 
Poles = 4.0 

Jm О о о at 
Km 09. 009007728 


!"Brushless DC Motor Equations" 


!"Rotor Ref Frame Voltages" 


2.*cost(thetar)*Vas/3. € 2.*cos(thera-c-2 "із УБИЕ 


Vqsr z 
+ 2.*cos(thetar + 2.*pi/3.)*Ves/3. 


2.*sin(thetar) *Vas/3. + 2.*sin(thetar-2.*pi/3.)*Vbs/3. & 


Vdsr = 
+ 2.*sin(thetar + 2.*р1/3.)*\/с<5/3. 


"Inverter Currents” 


cos(thetar)*iqsr + sin(thetar)*idsr 


ias 


cos (thetar - 2.*p1/3.)*igsr + sin(thetar - 2.*pi/3.)*idsr 


lbs 


ісз- -іа5-105 


Шисогог Бес сгаше Currents" 
ріаѕг - (-rs*iqsr - wr*Ls*idsr - wr*lamref + vaqsr)/Ls 
iqsr z INTEG (piqsr 77979) 

Bidsr = (-rs"idsr + wr’Ls*’igsr +» vdsr) /Ls 


А INTEG(Didsr | 0 0) 


'"Rotor Ref Frame Angle Thetar and wr" 
рыла << (Ше <- TL Jm 

wrm = INTEG(pwrm, 0.0) 

wr = (Poles/2.) * wrm 


pthetarm = wrm 


thetarm = INTEG(pthetarm : 0-0) 


!" Torque Zquations" 
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(Poles/2.) * ((3./2.) * lamref) * iqsr 


Te 


Tl = Km * (wrm*wrm) 
END ! * End of derivative’ 


DISCRETE reset 
thetarm = 0.0 


END 


DISCRETE reset2 
thetarm = thetarm + twopi 


END 


END ! “End of dynamic section’ 


END ! "End of program” 
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APPENDIX B. [ACSL CODE FOR NBC FAN MODEL WITH 120° INVERTER] 


PROGRAM NCB FAN MODEL (AAAV) 


INITIAL 


'"Preliminary stuff" 


MAXTERVAL maxt = 1.0e-7 !"max time step" 
MINTERVAL mint = 1.0e-7 !"min time step" 
CINTERVAL cint = 1.0e-5 !"communication interval" 
ALGORITHM ialg = 5S '"4th order Runge-Kutta" 
CONSTANT tstop = 0.007 !"adjustable stop: time" 
"---initial circuit constants---" 

CONSTANT rLz0.01 !"input line series К" 
CONSTANT LL-.000001 "input line series L” 
CONSTANT Cdc-.036 !"input capacitance" 
CONSTANT Vdc-28.0 !'"input voltage" 
!"---parameters for determining circuit status" 

CONSTANT hold = l.e-l i"holdrng current ror diodes: 
CONSTANT eps = 0.7 !"diode drop" 

!"Brushless DC Motor Constants" 

CONSTANT Ls - .000005 '"self inductance" 
CONSTANT rs = 0.003 "stator resistance" 
CONSTANT lamref = 0.006 . ! "magnet parameter" 
CONSTANT Poles = 4.0 ' | !"number of poles" 
CONSTANT Jm - 0.00011711 !"motor-load inertia" 
CONSTANT Km = 0.000001128 |! "TL2Km*wrmé&wrm* 


"---logical variables for transistor conduction status" 


ECGICAL TECN T2? о Td rS UPS 
"---logical variables for diode conduction status" 


LOGICAL Di 92-05% 70427705. DG 


"---logical variables derived £rom sensorless signals" 
LOGICAL Tla, T2, T32, Tea Toa, тоа 


"_---switch initial conditions----" 
"_---- -- all devices off" 
т = а се 

T2 = alse 

а = заїсе 

ma = false 

pce False 

Т6 = „false 

DI = .false 

D2 = .false 

D3 a se 

D4 false 

DS сасе 

D6 = .Zalse 
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"---all sensorless signals initially false" 


Tla = .false. 
T2a = .false. 


T3a = .false. 

T4a = .false. 

T5a = .false. 

T6a = .false. 

"---logical signals initiate integration of phase voltages" 


LOGICAL flaga, flagb, flagc 
flaga = .TRUE. 


flagb = .TRUE. 
flagc. s ОТАПЕ. 
"---logical symbol to switch over to sensorless signals" 


LOGICAL NOSENSE 
NOSENSE = .TRUE. 


END ! “СЕ паста)” 


DYNAMIC 


*.---simulation termination condition” 
TERMT (t .GE. tstop) 


DERIVATIVE 


'"constants needed for manipulating angles” 


CONSTANT pi=3.14159265 
twopi -"2 И ЕН | 
CONSTANT alpha = 1.047198 !"60 degrees" 


PROCEDURAL(thetar = thetarm) 
'"--~shift rotor electrical angle by alpha" 


thetarx= (Poles/2.) * tnetarm + alpha 


!"keep angle between 0 and twopi" 

IF (thetarx .GE. twopi) TEEN 
thetarx= thetarx - twopi 

ELSEIF (thetarx .LT. 0.0) THEN 
thetarx = thetarx + twopl 

ENDIF 

thetarx= MOD(thetarx, twopi) 

!"recover thetar" 

thetar = thetarx - alpha 

END !" End of Procedural" 


airing of SWicCches: 
1"1f NOSENSE 1s TRUE, use rotor position signals for 
І"1Е NOSENSZ 1s FALSE, use open-clrcult-derived signa 


TINO SENSE) К НЕМ 

Te Eosrarxc os ND DL AND (ene sos ЕЕ TEN 
Т1 = .єхмче. 
т< гат5ет 
T3 = .false. 
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BESETE 


ELSEIF 


ELSEIF 


ELSEIF 


EBSEIF 


ENDIS 
ELSE 


ОСКЕ false. 
Т5 = .false. 
TO = .true. 

( (chetarx .GE. pi/3.) АМО (thetarx ЕТ 2 руз у у ты 
TL = brune. 
T2 = .true. 
T3 = .false. 
T4 = . false. 
TS = .false. 
T6 = .false. 

((chetarx .GE- (2.*pi/3.]] -AND. (thetarx oto) teen 
Tl s .false. 
T2 - .true. 
T3 = .true. 
T4 = .false. 
TS = .false. 
T6 = .false. 

((thetarx .GE. pi) .AND. (thetarx .LT. (4.*pi/3.))) THEN 
Tl = .false. 
T2 = .false. 
T3 = true. 
T4 = „true. 
T5 = .false. 
T6 = .false. 

((thetarx .GE. (4.*p1/3.)) .AND. (thetarx .LT. (5.*pi/3.))) THEN 
Tl = .false. 
T2 = .false. 
ту = .false. 
T4 -= true. 
TS = .true. 
T6 = .false. 

((елеса СЕ (5910030) ) -ANDi thetarx т оо) PEST 
Tl = .false. 
Т2 = .false. 
T3 = .false. 
та = са се: 
TS stus. 
T6 z rue. 

Тї = тта 

2 =" 12а 

T3 = T3a 

T4 = Tda 

TS t=) oa 

75-0792 

!"---Determine diode conduction status" 

["“---establish the open-circuit voltage" 

Масос = (17722722) = ((3.* lamre=*wrecos (theter))/2.)) 

Vbgoe = (УЕ) 663. "Јатгегти сог(ецевах = (2. =>: и | 2) 
(\(Widey 2.) “+ ((3."*lemeeof"we cos (cnetar % (27522 3 0072000 


Vecoc 


t"diode logic for PHASE LEG A" 
1*3f either transistor is gated, then either the transistor 
" or diode will be conducting" 
IF ((T1) -OR BA) THEN 
D1 = .Еа15е.. 
D4 = .false. 


!"if the leg transistors are not gated:" 
!"if ias»hold, then the lower diode is on, upper is off" 
ELSEIF (ias .GT. hold) THEN 

D4 .true. 

D1 . false. 


!"if ias«-hold, then the upper diode is on, lower is off" 
ELSEIF (ias .LT. -hold) THEN 

D -erüe. 

D4 .false. 


"ТЕ open circuit voltages are between Vdc - eps and -eps" 
'"then both diodes are off" 
ELSEIF ((Vagoc .LT. (Vde + eps)) .OR. (Vagoc .GT. -eps)) THEN 
Dl = .false. 
D4 = .false. 


'"if open circuit voltage is above vdc + eps, upper diode is made" 
зүү On: 
ELSEIF (Vagoc .GT. (Vdc + eps)) THEN 

D1 = .true. 


!"otherwise the lower diode turns on" 


ELSE 
РА = Ене 


SENDIT 


'"diode logic for PHASE LEG В" 
ТЕ ((T I OR T6] T THEN 

Dio false. 

D6 = .false. 


ELSEIF (ibs .GT. hold) TREN 
D3 = .false. 
Dos Crue.: 


ELSEIF (ibs .LT. -hold) THEN 
Dis “смес 
D6 = .false. қ 


ELSEIF ((Vbgoc .LT. (Vdc * eps)) .OR. (Vbgoc .GT-. -eps)) THEN 
D3 = .ralse. 


be = 2156 


ELSEIF (Vbgoc .GT. (Vdc + eos)) THEN 
ра = Muere” 
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ELSE 
D6 = true. 


ENDIF 


|"алоде Јоде Еос PHASE БЕС С” 
IF (TS) eon (7279 THEN 

D5 = “false. 

D2 = .false. 


ELSEIF (ics .GT. hold) THEN 
DS .false. 
D2 . true. 


ELSEIF (ics .LT. -hold) THEN 
D2 - .false. 
DS = Crue: 


ELSEIF ((Vegoc .LT. (Vdc + eps)) .OR. (Vegoc .GT. -eps)) THEN 
D5 . false. 
D2 . false. 


ELSEIF (Vcgoc .GT. (Vdc + eps)) THEN 
DS- p erue. 


ELSE 
D2 - erue . 


ENDIF 


'"Procedural for finding stator voltages" 


AOS AL (Vas, Vos, Ves, Lili = Tl, T2, T3, T4, TS, TG, Di 52. 5. РА z 
DS, Bo, las, ibs, ics) 


! " Finding pole voltages and currents for A-pnase" 
'"upper device is on" 
ШЕН (ТЈ) .OR. (D1)) TREN 

Vag = Vdc 

iia = las 


!"lower device is on" 

ELSEIF ((T4) .OR. (D4)) THEN 
Мас = 0.0 
iia 0.0 


'"leg 1s open circuited" 


ELSE 
(ас л) + ((3.*larm ef*wr*cos(checaa a2 l) 


0.0 


Vag 


шла 


ENDIF 
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! " Finding pole voltages and currents for B-phase" 
EF (CTS ORC ЕТІНЕН 


Vog = vdc 
ilb = wos 
ELSEIF ((T6) (Өн (067) THEN 
ура = 0.0 
1152-200 
ELSE 
Vbg = ((Vde7v2.-) + 03. *lamref*wr*cosdftbetar исро 10 
ilb = 020 
ENDIF 


! " Finding pole voltages and currents for C-phase" 
ТЕ ((Т5) ОК, (D53) THEN 

Vcg = Vdc 

IIC ics 


ELSEIF ( (T2) F. OR D2) THEN 


Vcg = 0.0 
СЕ ото 

ELSE | 
Veg - (iVdc/2.) * ((3.*lamref£*wr*cos(thetar + (2. tp ор 
21e 050 

ENDIF 


I" Establish the phase voltages" 


vas = 0255Vag 30) 2 Vbg7 39-0 Veg 3) 
МОБИ ИСС tvboc4s)c-svag^s) -c-(Vco23:9) 
Vieew= milo. *“VE07o8) = (масл) = (VbgqgZ3.)) 


!" Find the inverter input current" 
iili = ila + ID + LITE 
!" End of Procedural" 


"evaluate the neutral voltage" 
vng - (vag + vbg + vcg)/3.0 


!"C-Phase Sensorless Signal derivation" 
CONSTANT we = 10000.0 !"low-pass filter frequency" 


CONSTANT tau2 = 0.004 
!"calculate the terminal voltage: should be Vcs" 


мајЕтс = уса = ла 4 


I" lowpass filter to eliminate abrupt transitions" 
CONSTANT xfIcic = 0.0 

рХЕІС = Ewex КоЛ We vdi iC 

хЕје = INTEG(pxrie T x Ie C) 
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!"activate integrator following first peak" 
SCHEDULE intone „XP. (16.0 - xflc) 
IF (flagc) then 

pxf2c = 0.0 


ELSE 
pxf2c 0552500: Doce 


ENDIF 
xf2c = INDTEG et 2e7 75070) 


apply result to third filter to eliminate dc componenti 


русх = (хЁ2с - wex)/tauz 
vex - ІМТЕСІруел 90) 
vco = -vcx * xf2c 


!"do same for the A-Phase" 
vdiffa = vag - vng 
CONSTANT xflaic = 0.0 
pxfla = -wc*xfla + wc*vdiffa 
xfla = INTEG(pxfla, xflaic) 
SCHEDULE intona .XP. (16.0 - xfla) 
IF (flaga) then 
pxf2a = 0.0 


ELSE 
29500 0556s da 


pxf2a 


xf2a = INTEG(pxf2a, 0.0) 
= (xf2a - vax) /tau2 
vax = INTEG(pvax, 0.0) 
= -vax + xf2a 


!"do-same for the B-Phase" 

ue ffb - vbg - vag ; 

CONSTANT xflbic - 0.0 

ЕХСІрБ = -weixflibrr woc*valrfb 

ЕТЮ = INTEGIpxEED, XETDbre) 
SCHEDULE intonb .XP. (16.0 - xflb) 
IF (flagb) then 


pxf2b = 0.0 
ELSE 
рхЕ25 Ss 7250020" xElb 
ENDIF 
БОБ - INTEG(pxr2b, 0.0) 
pvbx - (xf2b - vbx)/tau2 
mox s INTEG(pvbx, 0.0) 
уро = -vbx + x£2b 


vesoc = lamref*wr*cos(thetar + 2.*p1/3.) 


Link Dynamic Equations" 
!"removed for simulation development convenience" 


ош = (Vác - rL*IL - Vcap)/IL 

гір = INTEG(piL , 0.0) 

!'pVcap - (iL - iii)/Cdc 

BU cap = INTEGÍpVcap 7233.9) 

Усар - 28.0 "assume hard input voltace" 
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"Brushless DC Motor Equations" 


!"Rotor Ref Frame Voltages" 
Vgsr s 2.*"cos(tHeEar)*Vas/3. * 2.*cos(tbhetar-2 *53/3-) *U S MC 


+ 2.*cos(thetar + 2.*pi/3.)*Vcs/3. 


Vdsr = 2.*sin(thetar) *Vas/3. + 2.*sin(thetar-2.*pi/3.)*Vbs/3. & 
+ 2 WSim(thetar + 2. *pi/3.)*Ves/3. 


! "Phase Currents" 
ias = cos{thetar) *iqsr + sin(thetar) *idsr 
ibs = cos(thetar - 2.*pi/3.)*iqsr + sin(thetar - 2.*pi/3.)*idsr 


ics= -ıas-ibs 


'"Differential equation for the machine" 
piqsr =(-rs*iqsr - wr*Ls*idsr - wr*lamref + vqsr)/Ls 
así c PNPSSEOYGSrT 20950) 


pidsr = (-rs*idsr + wr*Ls*iqsr + vdsr)/Ls 
idsr = INITEC (pidsr . 0.0) 


I" Torque and mechanical dynamics" 
Te = (Poles 2.) waa 34/2.) cameo ase 
Tigo sm =) (Wr wim) 


pwrm - (Te - Tl)/ Jm 
wrm - INTEG(pwrm, 0.0) 
wr = (Poles/2.) * wrm 


pthetarm = wrm 
thetarm = INTEG(pthetarm , 0.0) 


!"use sensorless signals to derive switching signals" 
ТЕ (мес LT О с SAND. (уво -bT. 0.01 THEN 
Tla = .FALSE. 


T2a = .TRUE. 
T3a = .TRUE. 
T4a = .FALSE. 
[5a 2 .FALSE. 
T6a = .FALSE. 
ELSEIF ((vbo .GT. 0.0) .AND. (vao .GT. 0.0)) THEN 
Tia = .FALSE. 
T2a = .FALSE. 
Dia = TRUE. 
T4a = .TRUE. 
T5a = .FALSE. 
T6a = .FALSE. 
ELSEIF ( (vao ГТ „АМ О. (veo т. 0020 ТЕБЕМ 
Т1а = .FALSE. | 
T2a = .FALSE. 
Tia = .FALSE- 
Т4а = . TRUE. 
Sa "=e RUE... 
T6a = . FALSE 
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ELSEIF ((vco .GT. 0.0) .AND. (vbo .GT. 0.0)) THEN 


Tla = .FALSE. 
T2a = .FALSE. 
T3a = .FALSE. 
T4a = .FALSE. 
Тба = .ТАПЕ. 
Ta UTE 
ELSEIF ((vbo .LT. 0.0) .AND. (vao) LT 0 О THEN 
Т1а = DUE 
T2a = .FALSE. 
T3a = .FALSE. 
T4a = .FALSE. 
TSa = .FALSE. 
Тба =. TRUE. 
ELSETE ((мас Ст. 0.0) SAND: (vco .GT. 0.0)) THEN 
] Tla = .TRUE. 
T2a = .TRUE. 


T3a = .FALSE. 
T4a = .FALSE. 
TSa = .FALSE. 
T6a = .FALSE. 


ENDIF 


END ! " End of derivative” 


DISCRETE intona 
flaga = .FALSE. 
END ("of discrete inton" 
DISCRETE intonb i 
flagb = .FALSE. 
END ("OE discrete inton” 
DISCRETS intone 


flagc = .FALSE. 
END !"of£ discrete inton" 
! "End of dynamic section" 


"End or program" 
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APPENDIX C. [ACSL PLOTS FOR NBC FAN MODEL WITH 180" INVERTER] 


Machine parameters are those from Table 5.2. 
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APPENDIX D. [ACSL PLOTS FOR NBC FAN MODEL WITH 120" INVERTER] 


Machine parameters are those from Table 5.1. 
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APPENDIX E. [MATLAB CODE FOR ULTRA CAPACITOR MODEL | 


%%8%%%%%%%%%%%%%%%%%%%%9%%%%%%%%5%%%%%%%%%%%%%%%%%%%%%%%%%59%%%%%%%%%558%%%%9% 


&Capt Gabriel Beltran % 
%18 Әап 00 % 
$Ultra Capacitor Model % 


599499%99%%219%%99%994%9%%9%%99%%99%99%499%9199%9%%99%9%%99%9%%99%99%9%9%% 
$ Vendor Data 


figure(1) 


freq = |.01,.0159,.0251,.0398,.0631,.1,.1585,.2512,.3981,.6310,1. 0/1 58468 .. 


2,5119, 32 ӘБ 1; 6. 3096, 10.0,15.84839,25. 1189, 39.6107 76502) 
210005051557 4893,251.1887,:398.1072,630-95 24 000 О: 


mag = [23.4411,14.8375,9.6032,6.1331,3.93687 2.533997 664 7 PONE E ED 
0541805,093291,0.2325,0.1692,0.1275,0.0991,09079] 50 00399 050 P RE 
0: 0446,0.0384/0.0337,0.0300,0.0270; 0.0246,0.0226 070210): 


$ Plot 
semilogx(freq,mag) 


title('Impedance vs. Frequency!) 
xlabel (' Frequency (Hz)') 
ylabel (' Impedance (Ohms)') 


$ Transfer Function Simulation 


figure (2) 
а = 5.350; 


b-input('input b') 


[b, b*a]; 


num 
den = [1.0, 0.0]; 
magus - abs(freqs (num,den, freq) ); 


3 Plot 

semilogx (freq,magus) 

title ('Impedance vs. Frequency') 
xlabel('Frequency (Hz)') 
ylabel('Impedance (Ohms) ') 


% Overlay Plots 


figure (3) 

semilogx (freq,mag, '*',freq,magus,'.-') 
title('Impedance vs. Frequency') 
xlabel('Frequency (Hz)') 
ylabel('Impedance (Ohms)') 
EEexbt('Actual = tet) 

gtext('Model = .-.-.-') 
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APPENDIX F. [MANUFACTURER DATA SHEET FOR OPTIMA 800U] 


| Battery Mode: 8000 

Partnumber 8004-003 

Mo n inal Voltage: 12 volts 

NSN:6140- 01-3742243 

2” ph оп: pw power, sealed, lead acid, engine starting battery. 








|High риту jead-tin nally Wound cell configuration using utilizing proprietary 
ISPIRALCELL technology. 


Polypropylene 


Case: Dark Gray 
Cover: "Optima" Red 


D pt 4 
S НЕ, 
m ais E 
4 > M ` 
к; ES 
ж 2. ес 2 ` 
G 2$ ү: 
"i FAR =: y | 
| т 
14 PE y 
A Jf 
р = 
> 
~ 
ПШ itera 


Standard ¡Metnc 

10" 224 mm 

Dn. 6” 1722 mm 

! Een E. 78" 4981 mm(height at the foo of the terminals) 

| {тит ТТЕ СТТН 17.6 kg 

om pinal Configuration: SAE / BCI automotive and GM style side terminal (3/8 - 16 UNC - 2B, 
Я “ш 











: > n nic Voltage гү Charged): |12.8 volts 
(fema! al resistance (fully charged): |0032 ohms 


26 Ah (C/20) 50 Аһ (С/2) 





E BCI120 minutes 
Beo . .. | amp ciseharge, 00°F 267°C) o 0Svolecom| — 
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8% amps 


Lg 1000 amps 


| | 50 cycles 


7 i | E А Cy e (BC 10,000 - 12,000 cycles 


са 
Thek 


TC lowing charging ај are recommended to ensure a long battery life: — 
(АМЕ Ways Use a voltage regulated charger, with voltage limits set as described below.) 


ise batteries are designed for engine starting applications. They are not recommended or 


= = ЖЕЛП 


ed for use in deep cycle applications. 


Maximum Fer 19.0 volís (regul "T 
Maximum current: 60 amps. 
Maximum recharge time: 2 hours. 

All limits must be strictly adhered to 
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A f lis | ger safety glasses when working with batteries. 


Д ісі suse a voltage regulated battery charger with limits set to the above ratings. Overcharging 
салса use the safety valves to open and battery gasses fo escape, causing premature failure. These 
passe is are flammable! You cannot replace water in sealed batteries that have been overcharged. 
Any battery that becomes very hot while charging shouid be disconnected immediately. 


pret à 


lit ng and Transportation information: | 

Optima batteries can be shipped by AIR. The battery is non-spillable and is tested according to 
ICAO Technical instructions DOC. 9284-AN/905 to meet the requirements of Packing instructions 
| No: 06 andis classified as non-regulated by [ATA Special Provision A-48 and А-67 for UN2800. 
Те тігі must be protected from short circuit. 


Manu acturing Location: 

5 pima Batteries, inc. 

17500 East 22nd Avenue 

Auro fa, CO 80011 

| Int ed States of America 
Phone: 303-340-7400 

ED. 503. M 7474 


Ji: jarge versus temperature (400 Tos. 
Sel I Discharge: 
El i gineering Drawing: 


па Batteries, Inc. 
ici Specifications: Model $00U 
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